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We show that the temporal fluctuations AH(t) of the threshold driving field H(t), which triggers 
an avalanche in slowly driven disordered ferromagnets with many domains, exhibit long-range corre- 
lations in space and time. The probability distribution of the distance between successive avalanches 
as well as the distribution of trapping times of domain wall at a given point in space have fractal 
properties with the universal scaling exponents. We show how these correlations are related to the 
scaling behavior of Barkhausen avalanches occurring by magnetization reversal. We also suggest a 
transport equation which takes into account the observed noise correlations. 



I. INTRODUCTION 



Field driven disordered ferromagnets at low tempera- 
tures exhibit Barkhausen (BH) noise, a very important 
physical phenomenon, which is used for noninvasive char- 
acterization technique in commercial alloys. It has been 
recognized that measured BH noise exhibits scale free be- 
havior when only the external field is varied for all values 
of disorder and driving fields in real experiments (a 
short summary of the experimental data can be found in 
Rcf. |5j). This is in sharp contrast to some theoretical 
conclusions Q emphasizing fine tuning of the strength 
of disorder to a single critical value. It has been also 
attempted to understand the dynamics of domain walls 
which results in BH avalanches in terms of the dynamics 
of a sandpile model |?],|| (for a recent review of sandpilc 
models see ||) and of the models of interface depinning 
LO 12 1 . It should be stressed that in BH noise the 



presence of disorder plays an important role via pinning 
of domain walls. It remains to be understood how the 
motion of domain walls is affected by the pinning. 

Different scaling behaviors of BH noise can be at- 
tributed to the domain structure, which is related to an- 
nealing and type of impurities, and nucleation and coa- 
lescence of domains, as well sa varying driving conditions. 
Annealing the samples in the applied anisotropic stress or 
in the magnetic field leads to a characteristic structure of 
extended system-size domains with a 180° domain walls 
parallel to the anisotropy axis |T^,|lC],[llJ . Numerical sim- 
ulations with extended domain wall in two dimensions 
p2]| led to the conclusion that different scaling behavior 
can be expected in two limiting cases: (a) When disorder 
is weak BH response is dominated by motion of a single 
(extended) domain wall; (b) For strong disorder a mul- 
tidomain structure occurs with many competing domain 
walls. This picture is in a qualitative agreement with 
experimental results in stress annealed Fe-B-Si [|l0| and 
Fe-Co-B [[ll] alloys. Therefore the two different univer- 
sality classes can be related to surface and bulk critical- 
ity, respectively. Numerical simulations starting from an 
uniform ground state p| p^Jl5[ |, on the other hand, do not 
take into account extended domain walls, and thus cor- 



respond to the behavior at high-disorder. In both cases, 
however, the origin of scaling has not been fully under- 
stood. 

Recently the exact results of the random-field Ising 
model on the Bethe lattice Jl(| show that the avalanche 
distributions at a fixed driving field have finite cutoffs. 
However, an infinite cutoff appears for a range of disor- 
der strengths A < A c if a distribution is integrated over 
the hysteresis loop. The integration thus involves an in- 
finite jump in magnetization (when system size L —> oo) 
at a critical value of driving field. On the other hand, for 
strong disorder A above A c the avalanche distributions 
remain exponentially damped for all fields [jl6| . These re- 
sults encourage further study of the field-integrated dis- 
tributions to elucidate the role of driving field in the ap- 
pearance of the scale-free BH noise. 

In this work we study properties of the threshold driv- 
ing field which triggers Barkhausen avalanches in the 
multidomain structure, which is generated in the linear 
part of the ascending branch of hysteresis. When the 
system is slowly driven by increasing the external mag- 
netic field H with time, even by continuous field changes, 
there exists a threshold field that equals weakest pinning 
force of a domain wall in the system and thus starts an 
avalanche. In the numerical experiment we can adjust 
the driving field updates to the weakest local field in a 
system (so-called infinitely slow driving), and thus we can 
examine precisely the properties of the triggering noise. 
Surprisingly, we find that time series of such field updates 
are long-range correlated. Both Fourier spectrum of the 
threshold field fluctuations AH(t) = H(t + 1) - H (t) 
and distribution of distances of the successive avalanches 
which that field triggers appear to decay with a power 
law. The exponents are (weakly) universal in a range of 
values of disorder where BH avalanches have a power-law 
distributions. In addition, we find that the distribution 
of trapping times of domain wall at a given point in space 
exhibits scaling behavior. We show how the long-range 
noise correlations and fractal properties of the trapping 
time distribution are related to the observed scaling be- 
havior of BH avalanches. 

The present study is motivated by recently observed 
pattern formation and activity correlations along a 
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driven interface |18|-|20|] in 1 + 1 dimensions stacked by 
random defects. In contrast to these models, here we 
have a 2-dimensional system with many interacting in- 
terfaces, and two separate time scales (compared to ex- 
tremal driving): slow time scale of field updates, and 
avalanche propagation time scale between two field up- 
dates. Despite of these differences, in both cases the in- 
termittent avalanche-like dynamics occurs, which we be- 
lieve is essential for the observed scaling behavior. 



II. MODEL AND SIMULATIONS 

We consider a simple model with disorder represented 
by local random fields h x , which appears from an original 
disorder via coarse-graining 
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where x = (xu,xx) and J x ,x' = 1 is a constant interac- 
tion between nearest-neighbor spins S x = ±1. A Gaus- 
sian distribution of h x is assumed with zero mean and 
width /. (Other types of disorder have been also con- 
sidered ||,|lj,|l^|). A domain wall of reversed spins is 
created along < 11 > direction Ell on the square lattice 
rotated by 7r/4. This model is motivated by the extended 
domain walls in stress-annealed samples, as discussed in 
Ref. |12[] . Periodic boundaries are applied in the direc- 
tion of wall and an open boundary at the opposite side of 
the wall. The system is driven globally — updated value 
of the external field is applied to all spins in the system. 
The dynamics consists of spin flips when the local field 
exceeds zero (see below). 

As discussed in Ref. |l2| , the fact that motion of the 
(11) domain wall has no energy threshold at vanishing 
disorder has several advantages. In particular, the wall 
dcpinning occurs along a line H c (f) ~ nf in the ( f, H) 
plane, where k is a constant (see Fig. 2 in Ref. [fjjfl). 
Therefore for the domain wall along (11) direction an 
infinitesimally small field H is sufficient to move and de- 
pin the domain wall for small disorder, H c (f) — > for 
/ — > 0. This important property of the model allows 
us to use much smaller lattice sizes compared to earlier 
studies 22 (| , where a large system has to be used in 
order to find a random field large enough to surmount 
the energy barrier 2 J in the case of wall along (10) direc- 
tion, or the nucleation energy 4 J in the case of uniform 
ground states. On the other hand, by increasing disor- 
der the distance between pinning centers decreases, and 
thus a critical disorder /* exists at which depinning is no 
longer possible. Instead, it becomes energetically easier 
to nucleate new domains in the bulk. The exact value 
of /* depends on the type of distribution [jl^,[l6|, and 
is still not known exactly. In Ref. [Q it was found us- 
ing a finite-size scaling analysis of the averaged interface 
velocity and avalanche distributions that for the present 
model /* = 0.62 within numerical error bars. 



The scaling properties of BH avalanches in low dis- 
order regime were discussed in Ref. ]l^ |. Here we are 
interested in the high disorder region, where the (11) in- 
terface remains pinned at all fields. Nevertheless, the ini- 
tial state with the (11) interface ensures that the ratio of 
the characteristic lengths <C 1 holds for the applied 
disorder values (thus the avalanche size cut-off sq <C L 2 ). 
(The characteristic length Ei represents distance between 
strong pinning centers |22[ for given strength of disor- 
der.) On the contrary, when ~ 1, the (11) interface 
moves, meaning that the selected system size L is small 
for given strength of pinning, and thus system is not in 
the high disorder regime. Thus, having the condition 
^i/I < 1 satisfied for each applied value of disorder, we 
can concentrate on the effects of disorder fluctuations by 
applying many configurations of random fields at fixed 
/ and L values. We use up to L — 768 and up to 10 3 
disorder configurations (i.e., up to 36 x 10 6 spins). 

Since we are interested in the properties of the trig- 
gering noise, which are related to distribution of disor- 
der and the actual dynamics of the system, we avoid 
any "short-cuts" which can speed the algorithm [53. 
We apply the natural slow algorithm, which consists of 
the following steps: Random fields are stored at each 
site of an L x L lattice; The system is searched for the 
minimum local field h m i n — minf/i^), where h l ° c = 
Y] i J X , X 'S X ' + h x + H = hi r _ x + H , and then the driv- 
ing field is set to exactly h min + r] (we use r] = 10~ 10 ), 
which thus triggers an avalanche at that site; A list is 
made consisting of the sites which may flip at the next 
time step (neighbor sites to flipped spins); The spins on 
the list are examined for flip and a new list is made; The 
process is continued until no more spins can flip, then 
next minimum local field is found. It should be stressed 
that all spins in a single list (spin shell) are updated in 
parallel, i.e., in a single time step, in analogy to paral- 
lel update in cellular automata models. In this way we 
have well defined time scale of avalanche evolution (inter- 
nal time scale) as the number of steps that the updating 
procedure goes before the avalanche stops. On the other 
hand, the slow (external) time scale is set by the number 
of driving field updates. 

It was shown in Ref. [j^] that for low disorder / < 
/* = 0.62 the extended domain wall moves through the 
system and depinns at a critical field H c (f), however, 
above the critical disorder /* the built-in domain wall 
remains pinned and many new domains of reversed spins 
are nucleated inside the system. In this region, corre- 
sponding to high disorder (or low tensile stress), a typ- 
ical structure of clusters which occurs in the linear part 
of the hysteresis is shown in Fig. 1. The theoretical value 
/* can, in principle, be related to a critical value of ten- 
sile stress a* ~ below which the domain structure 
changes, as also observed in the experiment in Ref. [fl0|| . 
More precise experimental data would be necessary in 
order to determine the exponent fi. A rough estimate is 
that fj, is given by the correlation length exponent at the 
transition /*, (i « v ~ 2.3 p2j. 
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III. FRACTAL PROPERTIES OF TRIGGERING 
NOISE AND DOMAIN- WALL TRAPPING 

Time series of the magnetization jumps corresponding 
to the individual avalanches are known as BH noise. In 
the high disorder region either a new domain is nucleated 
and grown, or already existing domain extended, a BH 
pulse is associated with motion of a domain wall from 
one position to a new one. The area between two con- 
secutive domain wall positions corresponds to the size of 
BH avalanche (measured as the area covered by a sin- 
gle BH pulse). In principle, a domain growth process is 
not "linear" but fractal, leading to the dynamic exponent 
z = 1.23, and fractal dimension of avalanche D = 1.88 
Jl2[ |. As usual, the dynamic exponent z and the fractal 
dimension D are defined via the scaling of characteris- 
tic duration and size of avalanches with the change of 
length scale, (t)x ~ L z , and (s)l ~ L D , respectively. 
Scaling properties are studied in detail (see set of scal- 
ing exponents in Ref. Jl2|). In particular, in the criti- 
cal region above /* the distribution of size of avalanches 
D(s) ~ s~ Ts and duration P(t) ~ t~ Tt scale with the 
exponents t s = 1.30 and r t = 1.47, respectively. The 
avalanche distributions integrated over hysteresis loo p fo r 
/ > /* can be scaled according to the scaling form [[L2| 

P(tJ,L)=t- T 'V(t(dfr,tL- z ) , (2) 

where 5f = (///* — 1) and (8f)~ v measures the cor- 
relation length due to the critical point at /*. For in- 
stance, when the condition L ^> (5f)~ u is satisfied, the 
avalanche cut-offs can be scaled with disorder strength 
(in the entire region of power-law behavior), leading to 
/* = 0.62 and v = 2.3 within statistical error bars (see 
fl2"f for details). It should be stressed that in this re- 
gion of disorder the spatial extension of an avalanche is 
determined not only by strength of pinning but also by 
the blocking by previous clusters (see Fig. 1). The block- 
ing effects dynamically alter the strength of pinning and 
thus influence the scaling properties. This illustrates a 
complex interplay between the disorder and dynamics, in 
contrast to, for instance, the equilibrium clusters in the 
random-field Ising model. In order to get an insight into 
these complex phenomena, we study next the distribu- 
tions of wall trapping times and distances between the 
points with weakest pinning, from which an avalanche is 
released. We find that both of these distributions exhibit 
a scaling behavior. 

The distribution of distances between initial points of 
the successive avalanches is shown in Fig. 2, where we 
distinguish between distances measured in the parallel 
(x||) and transverse(j_) direction relative to the direction 
of the wall. In fact, due to anisotropy of domain wall 
motion in this case, these distributions show different 
scaling exponents t\\ = 1.04 ± 0.03 and t± = 0.60 ± 0.03 
according to 

G(x h x±) ~xpg(x ± /xf) , (3) 



where (a — T \\/ T ±- Notice that for small distances corre- 
lations in both directions are almost equivalent. How- 
ever, a crossover to a distinct power-law behavior for 
G(x±) occurs at a distance (x m 30 in Fig. 2) which is 
presumably related to the characteristic size of avalanche 
at given disorder. The open boundary in the direction 
opposite to wall leads to the cut-off at large x± . It should 
be stressed that spatial correlations are sensitive to dis- 
tribution of disorder. Therefore, here we used a large 
number of samples in order to minimize scatter of the 
data due to disorder fluctuations and to clearly distin- 
guish between correlations in parallel and perpendicular 
directions. 

In the inset to Fig. 2 we show the Fourier spectrum 
of the threshold field fluctuations AH(t) measured on 
the external time scale. It exhibits two correlated re- 
gions corresponding to earlier and later times, respec- 
tively. Steeper slope in the inset to Fig. 2 vary from 
0.6-1.03, and flatter slope from 0.25 to 0.36, depending 
on /. The long-range correlations exhibited in Fig. 1 
show that the system evolves in such way that its next 
relaxation event depends on the history of the present 
state of the system. It is interesting to note that a 
similar time series of the fluctuation of magnetization 
AM(t) = M(t + 1) - M(t) represents the Barkhausen 
noise itself (looked at the external time scale) . An exam- 
ple of the BH noise signal is shown in Fig. 3. Note that 
on the external time scale duration of each elementary 
signal is equal to one, whereas the height of each ele- 
mentary signal represents corresponding avalanche size. 
Sizes of the successive avalanches are only weakly corre- 
lated in time. However, a time derivative of the signal, 
namely a(t) = d[AM(t)]/dt, representing acceleration of 
a domain wall at each field update, shows certain cor- 
relation properties. In the inset to Fig. 3 we show the 
Fourier spectrum of the numerical derivative of the the 
signal, which is shown in main Fig. 3. Average slope of 
the curve is close to one. 

It should be noted that the average number of 
avalanches occurring at fixed strength of disorder in- 
creases with system size as (n a ) ~ L 2 / ' (s) ~ 
L^- D ^ T ^<j}(L(Sf) u ), where (j}(L(Sf) u ) is unknown 
scaling function. This implies that in the theoretical case 
of infinitely slow driving the average size of the field jump 
(AH(t)) decreases with L as (Aff(t)) ~ H sat (L)/{n a ) ~ 
L-°- 6S H sat (L)/(t)(L(SfY), where H sat is the saturation 
field. In real experimets the size of field jumps are re- 
stricted by the driving frequency, however, number of 
avalanches detected varies with the size of pick-up coil. 
This problem requires more detailed theoretical and ex- 
perimental investigation p5| . 

Further understanding of the role of disorder in the 
dynamics can be achieved by considering the time inter- 
vals that a domain wall resides at a given point in space 
(trapping time T trap ). We calculate the distribution of 
trapping times T trap of a domain wall at a given site, 
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which is determined as the time interval since the do- 
main wall is pinned at a site (x\\,x±) until it eventually 
moves away from that site. In this way, T tra p measures 
the time interval between two successive activities at that 
site. Notice that T trap in this case is somewhat different 
from so called first-return time in 1-dimcnsional interface 
p8| . Here we have 2-dimensional system with many in- 
teracting interfaces. Another similar example is trapping 
of grains in rice-pile models |^6| . A reasonable time scale 
to measure T trap is the number of field updates (external 
time scale) , since in the zero-temperature dynamics dur- 
ing an avalanche a spin at a given site is cither fixed or 
flips only once. 

The trapping-time distribution is shown in the inset 
to Fig. 4 for various values of disorder / > /*. Between 
the lower cut-off T (below which all trapping times are 
equally probable) , and an upper cut-off due to lattice size 
L, the trapping time distribution shows a short region 
with correlations. In particular, for large trapping times 
we can determine the slope as r trap = 2.3 ± 0.1. One can 
understand that a metastable configuration of domain 
walls corresponds to the system residing in a "pocket" 
of the fractal free energy landscape. A domain wall may 
move from a site (x|| , x±) as a part of an avalanche which 
started in the neighborhood of that point, corresponding 
to a local reorganization of the landscape near a shal- 
low minimum (for a given value of the external field). 
When the system resides in a deeper minimum, however, 
it waits for larger driving fields (i.e., trapping time in- 
creases) or for a more global reconstruction of the land- 
scape, which occurs with smaller probability compared to 
the one discussed above. This leads to the large slope of 
the distribution P(T tra p) ~ T tr ^ p ap at large T trap . It is 
interesting to note that a similar slope for large trapping 
times was found in the case of rice-pile model , where 
trapping of grains are considered, and it was argued that 
the slope is related to the roughness exponent of the rice- 
pile surface (see detailed analysis in Ref . [^6| ) . Analogies 
between an interface motion and rice-pile model have 
been established in the literature (27). Considering the 
avalanche exponents the analogy also applies for the BH 
noise with an extended domain wall in the limit of low 
disorder E^]. Notice that, in contrast to ricepile model, 
the range of power-law behavior in BH noise depends 
not only on the system size, but also on the strength of 
disorder, which restricts spatial extension of avalanches 
to s < s max <C L 2 and their durations to T < T max . 
Therefore, in this region of disorder we have rather small 
range of the power-law behavior. However, the scale- 
free behavior of the distribution of trapping times can 
be demonstrated via finite size scaling analysis when dif- 
ferent system sizes are used. We find that the following 
scaling form 

P(Tt rap , L) = L a V{T trap L~ ZT ) , (4) 

applies with the exponents which are weakly dependent 
on disorder a = — 0.4±0.05 and zt = 1.66±0.05 (see Fig. 



4) . By increasing the system size L the trapping times in- 
crease, leading to larger cut-offs of the distribution. We 
notice that the lower cut-offs To also increase, leaving 
a limited correlation range. However, the cut-offs scale 
nicely with T, as shown in Fig. 4. In the inset to Fig. 4 
we show how the trapping times distribution varies with 
disorder. By increasing disorder the cut-off To moves to- 
wards larger values and the correlation region shrinks, 
corresponding to a lesser correlations in the dynamics, 
which also manifests in the decrease of the characteristic 
avalanche size. In the limit of infinitely strong disorder 
the dynamics becomes completely random, consisting of 
individual spin flips which align along a local random 
field. 

It is interesting to note that unlike the avalanche dis- 
tributions in Eq. (|^), the distribution of trapping times 
in Eq. (|J) as well as G(x) (|j|) do not show explicit dis- 
order dependence apart from a weak dependence in the 
exponents. In fact, the disorder effects in this case are 
included in the gradient of driving force AT (i.e., the in- 
crements of the external field AH, which are adjusted to 
the minimum local fields). For example, for the distribu- 
tion of trapping times we have in general 

P(T trap , L, AT) = L a V{T trap {AF) ZT / XF ,T trap L- ZT ), 

(5) 

where the increase in the driving force AT = H(t + 
Ttrap)— H(t) contains all the contributions due to interac- 
tion and random pinning which occurred during the time 
interval (t,t + T trap ), i.e., AT = £*+ T *™» h ir , x , which 
steam from the different sites in the system. According 
to the above results in Fig. 2, sum of these contributions 
has no characteristic scale. Therefore we may conclude 
that (AT) -1 / Af — » oo. Hence the first argument in the 
right hand side of Eq. (||) can be neglected compared to 
Ttrap/ L Ztrap , leading to the scaling form (|J), which is in 
agreement with the scaling plot in Fig. 4. 

Another interesting observation regards the dynamic 
exponent zt, which scales the tail of the trapping time 
distribution with the system size L. It may be related 
to the scaling of the length of the optimal path in 
strongly disordered medium, zt = Dop- In fact, the op- 
timal path between two points can be constructed from 
the most persistent sections of the domain walls in the dy- 
namics of BH avalanches. The length of the optimal path 
scales with the linear distance between the end points 
with the exponent Dop = top = 1-66 [p8p9| l. 

IV. TRANSPORT EQUATION WITH NOISE 
CORRELATIONS 

The long-range noise correlations are shown to be rel- 
evant for the scaling properties of the interface motion 
ppf . In the literature the noise correlations in the in- 
terface depinning are viewed as originating from an- 
other (external) dynamic processes [p0[ . In the case of 
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Barkhausen avalanches we see that the noise correlations 
appear as an intrinsic property of the dynamics when 
the system is driven infinitely slowly. The active section 
of domain wall can be represented by a surface h(x, r) 
which is pinned by quenched defects. The transport at 
the site (x, r) is described by the equation 



dh/dr = v\\dnh - 



v±&h- 



r)(x,t,h) 



(6) 



Here we distinguish total elapsed time r, which is de- 
fined as the sum of evolution times of all individual 
avalanches, on one hand, and t, which represents time 
scale of field updates, on the other. The quenched noise 
rj{x, t, h) is generated by varying the driving field in the 
steps which are adjusted to the minimum local strength 
of pinning, as discussed above. We expect that a dom- 
inate /i-dependence of the noise 77 can be expressed as 
r\ ~ p(x)dh, where p{x) represents the (anisotropic) lo- 
cal velocity of domain wall per field rate. In addition, 
time variation of the noise are related to the updates of 
the external field, and thus rj varies on the external time 
scale only. Therefore, we can write 



n(x, t, h) « p(x)(a + AH{t))(dh) , 



(7) 



where ao is a constant which depends on initial con- 
figuration and the hysteresis loop properties. The lo- 
cal interface velocity per field rate, p(x), is a random 
variable which is determined by the spatial distribution 
and strength of pinning associated with a given value of 
the driving field H, i.e., history of the system. There- 
fore p(x) is governed by a probability distribution, which 
can be deduced from the properties of the probability 
that an avalanche starts at distance x from the preced- 
ing avalanche, i.e., 



< p(x)p(x ) >= jG(x — x ) 



(8) 



It is assumed that the same type of correlations apply 
for the successive activities during the evolution of an 
avalanche, while the external field is constant. In the 
previous section we found that the probability G(x) ex- 
hibits long-range correlations in the perpendicular and 
parallel direction as G{x\\, x±) ~ xj^x^ ' 60 . We also no- 
tice that AH(t) ~ t~ e , where 8 is finite in the case of 
infinitely slow driving discussed above, whereas, 9 = in 
the case of finite driving rate in small steps pl| , where 
AH(t) = const. Notice that Eq. @ together with Eqs. 
(0) and (||) leads to an effective nonlinear term of the 
form jG(xh, x±)(ao + t~ e ) 2 (dh) 2 , which differs from the 
Kardar-Parisi-Zhang nonlinearity |3l| | by the power-law 
correlations in the coefficient. As discussed above, the 
origin of these correlations lies in the the dynamically 
varying pinning and blocking effects when a multidomain 
structure is slowly driven through the hysteresis loop. 

The relevance of the noise correlations for the uni- 
versality class of the interface depinning has been dis- 
cussed in the literature using dynamic renormalization 
group (RG) JjOl. Another example where correlations 



of random noise play an important role is represented 
by the scaling properties of river networks |2£||32|]. In 
that case the correlations of the type G ~ xjT 2+<5 x^_ 
were assumed, where 5 is an expansion parameter and 
the anisotropy exponent £ has to be determined self- 
consistently at a fixed point of the RG p8j. In the result, 
the scaling exponents do depend on the range of corre- 
lations 5. Similarly, we can expect that the long-range 
correlations of the form given in Eqs. (^) and (ph can 
be related to the universal scaling exponents of the BH 
noise. The exponents z, (q, which are defined above, 
and the roughness exponent x, which governs behavior 
of the dynamic variable h with the change of scale, can 
in principle be obtained by the dynamic RG applied to 
Eq. (||) with noise properties in Eqs. (7[^). Then the 
avalanche exponents r f and r s can be deduced using scal- 
ing relations in the critical region. The RG analysis of 
the transport equation (Q) with noise properties specified 
by Eqs. requires additional work and is left out of 

the scope of this paper. Here we discuss only the scaling 
relations between the exponents at RG fixed point, e.g., 
the dynamic exponent z, and the avalanche exponents. 

It should be stressed that such scaling relations are not 
universal and that they depend on the nature of the dy- 
namic process (see for instance ^,[32| for the case of river 
networks) . We argue that the scaling relation derived be- 
low are valid for BH avalanches in the high-disorder (or 
multidomain) region. The rational behind these scaling 
relations is found in the directed nature of the avalanche 
propagation. The evolution of an avalanche (see Fig. 
1) can be visualized as growth of a directed percolation 
(DP) cluster projected p3| to 1 + 1 dimensions, with extra 
dimension representing time axis. The fractal dimension 
of the equivalent DP cluster measured with respect to the 
time axis is then Du = D/z, where D and z are fractal 
dimension of BH avalanche and dynamic exponent, re- 
spectively, as defined above. Then D11 = 1 + (dp — Sdp, 
where (dp = z/2 is the anisotropy exponent measured 
with respect to time axis, and Sdp is the survival proba- 
bility exponent of the equivalent DP clusters. Notice also 
that for the directed dynamic processes r t = Du and thus 
t s = 2 — 1 /rj lU , which completes the set of avalanche 
exponents. In addition, the roughness exponent x can be 
related to the trapping time distribution as T tra p = 2 + x, 
as noticed in Refs. [^6|,^7). For instance, using the nu- 
merical values for D and z from Ref. JL2), the following 
values of the avalanche exponents are predicted by the 
above scaling relations: r t = 1.52, t s — 1.34, compared 
with 1.47 and 1.30, and x — 0.23, obtained by direct 
numerical simulations. 



V. CONCLUSIONS AND DISCUSSIONS 

We have demonstrated numerically the existence of 
long-range correlations in triggering noise which is in- 
trinsic to the domain wall dynamics in slowly driven dis- 
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ordered ferromagnets. Although the distribution of dis- 
order is initially uncorrelated, the spatio-temporal cor- 
relations develop in time. The appearance of these cor- 
relations can be attributed to the applied global driving 
in which always next weakest pinning force in the sys- 
tem is selected, on one hand, and to a finite extension of 
avalanches occurring between two consecutive field up- 
dates, on the other. Hence, attempts to reduce the do- 
main wall dynamics to the Glauber spin flips in the pres- 
ence of local random fields appears to be inadequate for 
the range of disorder where the cooperative avalanche dy- 
namics occurs. We suggest an alternative transport equa- 
tion which incorporates the observed noise correlations. 
The long-range correlations of triggering noise can be 
related to the universal scaling behavior of Barkhausen 
avalanches and to the transport properties via the frac- 
tal distribution of the domain wall trapping times. We 
believe that an analysis of the transport equation by 
the dynamic RG will contribute to understanding of the 
universality classes of Barkhausen avalanches by the in- 
finitely slow driving. By varying the driving conditions, 
however, these correlations are changed. This might be 
the origin of different scaling behavior of BH avalanches 
at finite driving rates, as observed in experiments [ pT| . 
Our results also suggest that technical procedures which 
speed the numeric algorithms in driven disordered sys- 
tems should be taken with care. In particular, an algo- 
rithm which alters the properties of triggering noise may 
lead to a different scaling behavior, which is unrelated 
with the original problem. 
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FIG. 1. Extended domain wall along the bottom line 
(shown in red), and unflipped spins (full blue color) and clus- 
ter of flipped spins developed at different stages of evolution 
(shown using a continuous color map [23]). Only 253 recent 
clusters are shown, with all older clusters shown in yellow 
color. L = 200 and /=1.1 . 
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FIG. 2. The probability that the subsequent avalanche 
starts at the distance x^,x± (bottom, top) from the point of 
the preceding avalanche for L = 192 and / = 1.0. The data 
are averaged over 1000 samples. Inset: Fourier spectrum of 
the driving field fluctuations AH(t) taken in one ascending 
branch of hysteresis for / = 1.0 and L = 192, vs. undimen- 
sional frequency. 
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FIG. 3. Sequence of Barkhausen pusles AM (measured in 
number of spins) recorded on the time scale of field updates. 
Inset: Fourier spectrum of the numerical derivative of the 
sequence in the main figure plotted against undimensional 
frequency. Fitted slope is 0.95 ± 0.03. 




FIG. 4. Inset: The distribution of trapping times of a do- 
main wall vs. trapping time Tt rap measured by number of 
field updates. The distributions for various values of disorder 
/=0.9, 1.0, 1.1, and 1.2 (left to right) are obtained at lattice 
768 x 768 averaged over 10 samples. The data are logarith- 
mically binned. Main: Finite size scaling plot of the trapping 
time distribution for /=1.0 and the linear lattice size L=192, 
384, and 768. The data for L = 192 are averaged over 1000 
samples, the rest of the data over 10 samples, and are loga- 
rithmically binned. 
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